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Abstract —In this paper, we formulate the interference align¬ 
ment (lA) problem for a multiuser multiple-input multiple-output 
(MIMO) system in the presence of an eavesdropper as a rank 
constrained rank minimization (RCRM) problem. The aim of 
the proposed rank minimization lA schemes is to find the 
precoding and receiver subspace matrices to align interference 
and wiretapped signals into the lowest dimension subspaces while 
keeping the desired signal subspace spanning full available spatial 
dimensions. To deal with the nonconvexity of the rank function, 
we present two convex relaxations of the RCRM problem, namely 
nuclear norm (NN) and reweighted nuclear norm (RNN), and 
transform the rank constraints to equivalent and tractable ones. 
We then derive a coordinate decent approach to obtain the 
solutions for lA schemes. The simulation results show that our 
proposed lA designs outperform the conventional lA design in 
terms of average secrecy sum rate. On the other hand, our 
proposed designs perform the same or better than other secure 
lA schemes which account for low interference and wiretapped 
signal power rather than for low dimensions of Interference and 
wiretapped signal matrices in the systems which achieve the 
perfect lA. 

I. Introduction 

Physical layer security (PLS) has recently attracted con¬ 
siderable attention due to its potential to improve the secu¬ 
rity of communication systems [1]. Different from traditional 
cryptographic methods which are implemented in the net¬ 
work layer, PLS exploits physical characteristics of wireless 
channels to provide secrecy. Secure PLS approaches have 
been developed for various communication scenarios, such 
as multiuser multiple-input multiple-output (MIMO) systems 
[2]-[4], relay networks, cognitive radio systems and other 
networks (see [5]-[7] and references therein). The transceiver 
designs to maximize the secrecy rate in the multiuser MIMO 
interference channels are of great concern. However, such 
design problems appear mathematically intractable due to their 
high nonlinearity and nonconvexity [2], [3]. Alternatively, 
interference alignment (lA) is one of potential techniques 
to increase the secrecy sum rate (SSR) in multiuser MIMO 
systems [8]-[10]. It is proven in [8] that it is possible for 
each user in networks to achieve a nonzero secure degrees of 
freedom (DoF) when using an lA scheme to design precoding 
matrices at each transmitter (Tx). The secure transmission has 
been shown to be feasible where the number of antennas at 
legitimate Tx and receiver (Rx) is greater than those of the 
eavesdropper [9]. The key idea of a secure lA design is to 


keep the desired signal free from interference and to offer 
a secure from the eavesdropper. To obtain interference-free 
desired signal, interferences at each receiver are aligned into a 
reduced-dimensional receive subspace, and to prevent desired 
signals from the eavesdropper, the wiretapped signals at the 
eavesdropper can be aligned into a proper subspace where their 
powers are minimized, or can be forced to zero [10]. 

The underlying secure lA problems are NP-hard and in¬ 
tractable. To handle these mathematical challenges, reference 
[10] introduced secure lA approaches which aim at minimizing 
interference and wiretapped signal power rather than reducing 
the dimension of interference and wiretapped signal subspaces. 
Recently, for the multiuser MIMO systems not in the security 
context, instead of minimizing the power of the interference 
signal, references [11], [12] proposed rank constrained rank 
minimization (RCRM) problems which ensure that the inter¬ 
ferences fall into a low-dimensional subspace and the desired 
signal spaces span all available spatial dimensions. Motivated 
by these works, for the multiuser MIMO systems in the 
presence of an eavesdropper, we adopt the RCRM framework 
in which the ranks of the wiretapped signal matrices and 
the subspace spanned by interference signals are minimized 
subject to full-rank desired signal space. In order to solve the 
NP-hard nonconvex rank minimization problem, we introduce 
nuclear norm (NN) and reweighted nuclear norm (RNN) as 
two convex approximations of the rank function, then propose 
heuristic lA algorithms to obtain near optimal solutions. Our 
experimental evaluations reveal that the proposed lA designs 
outperform the conventional IA design [13], and provide the 
same or better SSR performance than secure lA designs in 
[10] for the system in which the perfect lA can be achieved. 

Notations: Bold lowercase and uppercase letters represent 
vectors and matrices respectively. denotes the conjugate 
transposition of matrix X. X ^ 0 represents the Hermitian 
positive semi-definite matrix X. Id and 0 are respectively an 
identity matrix with d dimensions and a zero matrix with the 
appropriate dimensions. tr(.), rank(.), E(.) and [j:]+ are the 
trace, rank, expectation and max(x, 0) operators, respectively. 
||X||f is the Frobenius norm. A complex Gaussian random 
vector with means x and covariance is represented by 
X ~ CN{x,Rx)- amin{X) and ai{X) are respectively the 
minimum and the i-th largest singular value of matrix X while 

rank{X} 

ll-^ll* = S (T 2 (X) is the nuclear norm of matrix X. 

2=1 


II. System Model and Problem Formulation 



“ • • —► Wiretapped link 

Fig. 1. A system model of secure multiuser MIMO communication network. 


Consider a MIMO interference channel with K Tx-Rx 
pairs and an eavesdropper as depicted in Fig. 1. Without loss 
of generality, assume that the eavesdropper is the {K + l)-th 
Rx with the Nr^ antennas, each Txs or Rxs is equipped with 
the same number of antennas Nt and Nr, respectively. The 
fc-th Tx sends d data streams to the fc-th Rx. Such a system 
is denoted by an {Nt x Nr, Nr^, d)^ system [10]. Assuming 
that Hk,i € is the static flat-fading MIMO channels 

matrix from the £-th Tx to the /c-th Rx, the received signal 
yk S at the fc-th Rx for fc € /C = {1, ...,K -f 1} is 

given by 


K 

yk = '^Hk,eFese + nk (1) 

e=i 

where Fk € is the precoding matrix applied on the sig¬ 
nal vector Sfc G and Uk CA/" (O, is a complex 

Gaussian noise vector. The signal vector Sk is independently 
identically distributed (i.i.d.) such that EjsfcS^} = Id- We 
consider the scenario that a communication group involves the 
eavesdropper (the (iT-l-l)-th Rx) but the information is desired 
to be secret to him [10]. Thus, it can be assumed that the global 
channel state information (CSI) is available. 

The problem of our lA design for secure multiuser MIMO 
systems is to find the precoder matrices and the 

receiving signal subspace matrix which satisfy the 


three following conditions [10] 



Tank{WkHk,kFk) = 

d 

(2) 

W^Hk,iF, = 

0 ; ye^k,iGic 

(3) 

HK+i,tFi = 

0 ; G 1C, 

(4) 


where (2) enforces the useful signal to span all d dimensions, 
(3) and (4) ensure all the interference subspaces and wire¬ 
tapped signal subspace to have zero dimensions. 

Remark 1. The feasibility of a set of three linear lA equations 
above has been studied [10]. It is shown in [10] that perfect IA 
for an [Nt x Nr, Nr^,d)^ system is almost surely achieved 
under properness condition Nt — d > Nf, and Nr > Kd. If 
this condition is not satisfied, the system is called improper. 


For ease of presentation, we define the desired signal matrix 
Sk G interference matrix Jk G and the 

wiretapped signal matrix Se G ^ g 

respectively as follows 


Sk {Wk,Fk) =Wj!Hk,kF 


Jk {Wk,{Fe}l,^,^k) = 

Se = [{HK+uFe}t^ 


(5) 

(6) 
(7) 


Accordingly, (2), (3) and (4) can be restated in terms of ranks 


rank(5fc) = d, (8) 

rank(Jfe) = 0, (9) 

rank(5e) = 0, (10) 


for all fc G /C. To aim at obtaining d interference-free 
dimensions for every Rx and zero wiretapped-signal-space 
dimensions, we seek the proper precoding and receiving sub¬ 
space matrices to minimize the rank of interference matrices 
in (9) and the rank of wiretapped-signal matrix in (10) subject 
to the full rank of the desired signal matrix. Thus, the design 
problem of interest can be mathematically expressed as 

K 

min rank (Jk) + rank (Se) (Ha) 

}fc=l}fc=l k—l 

s.t. rank(5fc)=ci, € 1C. (Hb) 


Both the cost function and constraints in (11) are non-convex 
and, thus, problem (11) is mathematically intractable. There¬ 
fore, it is highly desired to find the efficient iterative algorithm 
to obtain suboptimal solution to (11) instead of the globally 
optimal solution. Specifically, we introduce convex surrogates 
for this cost function and the feasible solution set. This results 
another challenging problem is that, although the surrogate 
of the cost function is convex in either of the two sets of 
input matrices, the cost function is no longer convex if it 
is minimized over both sets at the same time. To handle 
this problem, we utilize a coordinate descent approach to 
alternatively minimize the cost function over the transmit 
and then over the receive subspace matrices iteratively. The 
constrain rank(5fc) = d can be replaced approximated with 
the following closed convex set, in the same manner as [11] 

SkhO, o-mm(S'fc) > e, ( 12 ) 


where e > 0. Note that, the orthogonality constraints on 
the precoding and receiver subspace matrices are omitted 
in problem (11). This is because that they can be linearly 
transformed to be orthogonal by using QR factorization, while 
the ranks of interference, desired signal and wiretapped signal 
matrices remain unchanged [11]. 


III. RCRM BASED lA APPROACH FOR SECURE MIMO 
Interference Channels 


A. NN approximation based lA algorithm 


To deal with the nonconvexity of the rank function, NN 
||X||* is used as a surrogate for rank(X) [14]. Problem (11) 
can be recast as follows 


min 

S.t. 


K 




* 


Sk ^ 0 ; tTmin{Sk) ^ 


(13a) 

(13b) 















Then, we solve the problem above by the coordinate descent 
approach as follows. 

Transmit precoder selection: When {W are fixed, 
for selecting {Fk]k=i^ (13) reduces to the optimization prob¬ 
lem Vp^ defined as 

K 

VT-. min ||5e|U 

s.t. Sk ^ 0 , crk) ^ r. 


Receive subspace selection: We use the solution 
of problem Vp^ as an input to select {Wk}^^i by solving 
problems defined as 


V, 


,NN 

Wu 


: min 
s.t. 


K 

k^l 

^ 0, (ymin (Sk) > e, 


The optimization problem above is convex and can be ef¬ 
ficiently solved using CVX toolbox [15]. We then feed the 
solution {Wk]k=i optimization problem back to Vp^. 

This process is continued over iterations. The step-be-step 
procedure is stated as Algorithm 1. 


Algorithm 1 : Secure NN lA Algorithm 


Inputs: d,Hk,e, V/c G ICU{K - 
where k is the iteration index; 

Initial variables: random matrix {F 


ipio) 


1}, \/i G 1C, K = 0, Kmax, 
k }f=i satisfied 


k /c 

while K < K 




do 


For fixed select by solving 


-ONN 


and orthogonalize 


For fixed {W), }k=i^ select {Fj. }k=i by solving 




NN 


and orthogonalize 


6 : K = K + T, 

1: Repeat steps 4-6 until convergence or when k reaches 

the maximum number of iteration Kmax- 

8 : end while 


The NN approximation in the cost function of (13) accounts 
for the sum of the singular values of the interference matrix and 
wiretapped-signal matrix, rather than the sum of the number of 
singular values. Thus, such an approximation may not result 
in the minimum rank of the cost function of (11). In order to 
make a tighter approximation of the rank function, we employ 
a different surrogate, namely RNN, for the cost function in the 
following subsection. 

B. RNN based lA algorithm 

The RNN approximation of a rank function was provided 
in the multiuser MIMO interference channels but not in 
the security context [12]. In this paper, we adopt the RNN 
approach to a secure multiuser MIMO system and develop 
an iterative RNN algorithm to reduce the wiretapped-signal 
and interference subspace dimensions. The basic idea of RNN 
approximation is using the function logdet(X + 61), known 


as a smooth surrogate of rank(X) where X is positive 
semidefinite matrix. The log-det type function is proved to 
possibly obtain low-rank solutions to linear matrix inequality 
problems for positive semidefinite matrices [14]. Using the 
RNN approximation, the RCRM problem (11) can be rewritten 
as 


min fl 

s.t. SkhO 

■) ^min {Sk) > e, 


(14a) 

(14b) 


where U = X) E log (cr*(Jfc) + 7 ) + E log(cr*(S'e) -f C). 

k—1 i—1 i—1 

de = mm{Ne, Kd), 7 and ( are the arbitrable small positive 
values to make the approximation resemble the rank function 
and ensure the stability for the log function. 


Since the cost function in (14) is concave, we apply a 
majorization-minimization (MM) algorithm to solve problem 
(14). In particular, at the K-th iteration, we find the upper bound 
of the cost function, then, minimize it to ensure that the cost 
function can reach the optimal minimizer over iterations. The 
upper bound of the cost function (14a) can be obtained by 
taking the first order Taylor expansions with respect to ai{Jk) 
and ai{Se) as follows 


log(CTj(Jfc) -f 7 ) < log +7) + 


CTjjJk) - f^ijj^k'^) 

(15) 


log (cr*(5e) + C) < log -f C) 


g,(5e) - cri{Si''^) 

(16) 


It should be noted that Se and their singular values are 
known at the K-th iteration. Hence, we iteratively minimize the 
upper bound given by 


K d 

EE 


k) 


E 




S.t. 

which is rewritten as 


mm 


K 


k^l 


S.t. 


^k E 0; k) E 


with condition Nr^ < Kd or 

s.t. Sh E ^m.in, iSk) > 


(^dexde jjjg vveight matrices which is defined as 



(17a) 

> e, 

(17b) 

Sell* 

(18a) 

; e. 

(18b) 


(19a) 

i e, 

(19b) 

and G 


in which G 


r(«) 


, A 


are the left singular vectors of 


X de 


are the left singular vectors of S, 


(k) 





when Nr^ < Kd ov the right singular vectors of 5^”^ when 
Nr^ > Kd, G is the diagonal matrix whose 

i-th diagonal element is equal to and 

G (j^dexde jj jjjg diagonal matrix whose i-th diagonal 
element is equal to = —- 7^5 -. The similar proofs of 

(TiiSe )+7 

these above transformations can be found in [ 12 ] and, thus, 
omitted. 

Now, we propose a RNN two-loop algorithm for the RCRM 
problem (14). At the K-th iteration, we create an inner loop to 
find optimal precoding and receive subspace matrices in (18) 
or (19) via the coordinated descent approach. The main steps 
in the m-th iteration inside the inner loop of our algorithm are 
stated in the following. 

Transmit precoder selection: By holding {Wk}^^i fixed, 
we select by solving the optimization problem 

defined as follows 

: min ^ 

S.t. Sk ^ 0, (TminiS ^ C. 
when Nr^ < Kd, or 

min ^ + ||5e$M||. 

s.t. Sk ^ 0 , fTminiSk) ^ r. 
when Nr^ > Kd. 

Receive subspace selection: We now set the solution 
of problem as an input to find receive 

subspace matrices {Wk}k=i solving the two following 
problem defined as 



K 

: min 

k^l 

s.t. 

Sk ^ 0'min{Sk^ ^ C, 


Each problem above is convex and can be solved by CVX 
toolbox [15]. This process continues until the cost function 
converges or m attains a specified maximum number of 
iteration rumax- We then update input and to the next 
minimization in (K+l)-th iteration. The step-by-step algorithm 
is described in Algorithm 2. 

IV. Simulation Results 

In this experimental evaluation, we run simulation for a 
(18 X 12,9, 3)^ and a (15 x 15, 9,3)^ system. Noise variances 
are normalized = 1. The Rayleigh fading channel 

coefficients are generated from the complex Gaussian distri¬ 
bution CA/"(0,1). We define signal-to-noise-ratio SNR = l^. In 
additional to the proposed lA designs, namely NN lA and RNN 
lA design, two other secure lA algorithms which minimize the 
power of interference and wiretapped signals are considered, 
such as the wiretapped signal leakage minimization (WSLM) 
[10] and zero-forcing wiretapped signal (ZFWS) [10]. It is 


Algorithm 2 ; Secure RNN lA Algorithm 
I: Inputs: d,Hk,i, Vfc G /C U {K -f 1}, Vf G /C, m = k = 0, 

^max? trijjiax^ 

Initial variables: 


9: 

10 : 

II: 

12 

13 

14 

15: 


random matrix 


p(o)H pW _ 

^ k ^ k — 


- / 

d -'d, —— -t, 


satisfied 


d and = Id', 


while K < Kmax do 
while m < mmax do 

For fixed {F^^^}k=p select by solv¬ 
ing and orthogonalize 

For fixed select by 

solving and orthogonalize ; 

Update m = m + T, 

Repeat steps 5-7 until convergence or when m 
reaches the maximum number of iteration TOmax; 

end while 

Output: and 

Evaluate ^nd from ( 6 ) and (7); 

,(k;+1 ) , /f,(K+l) 


Update E). 


and $ 


from ( 20 ); 


Update K = K + T, 

Repeat steps 4-13 until convergence or when k reaches 
the maximum number of iteration Kmax- 

end while 


noted in [ 10 ] that these two systems considered are proper for 
all these four lA methods, i.e.. Nr < 

Nr,. < Nt — d and Nr > Kd. We present the numerical 
results averaged over 200 channel realizations. We then plot 
and compare the achievable SSR of each lA design and the 
conventional lA design in which the security context was not 
taken into account. The channel capacity at the fc-th Rx, for 
k G JC, can be calculated directly from (1) as [16] 

TZk = log 2 \lNr + Hk.kFkF^ ( 21 ) 


K 

where ^ Hk/FiFf ^ + a^INr is the interfer- 

1=1,t^k 

ence plus noise correlation matrix in Eq. (1). The information 
leakage rate from the fc-th Tx to the eavesdropper can com¬ 
puted by 


FT = log 2 


'■Nr 


HK+i,kFkFk 


( 22 ) 


K 

where R^.k = Y. -ffA+i,r-Fr-Ff-fff+ 1 1 + (^llNr^ i® 

1=1,i^k 

interference plus noise correlation matrix at the eavesdropper. 
The fc-th Tx-Rx pair can obtain the secrecy rate given by 
Fs,k = \Fk — ■ Hence, the multiuser MIMO system 

can have the SSR 

K K 

Fs=J2 '^s,k = - nTT- (23) 

k=l k=l 

In the following experiments, we set the specific number of 
iterations for each lA design so that all lA algorithms take 
comparable times to run using MATFAB. For each simulation, 
we run 5 iterations of the NN lA algorithm, 3 iterations of the 
RNN lA algorithm, 50 iterations of the ZFWS lA algorithm 
and 25 iterations of the WFSM lA algorithm. To solve V2, 
where a = {F,Wk} and 77 = {NN,RNN}, we set e = 0.1 
and the maximum number of inner loops rrimax = 3. 








Fig. 2 plots the average SSR versus SNR for a (18 x 
12, 9,3)^ system. It can be seen that our proposed lA designs 
significantly improve the SSR as compared with the conven¬ 
tional lA design. Two proposed designs offer a slight SSR 
improvement when compared to the WSLM lA scheme while 
they provide the same SSR performance as ZFWS lA scheme. 
In this system, our lA algorithms appear to achieve the perfect 
lA and, thus, the RNN lA approach seems not to show the SSR 
enhancement as compared with the NN lA approach. 

The average SSR for a (15 x 15,9,3)^ system is plotted 
in Fig. 3. It can be observed that our proposed lA algorithms 
offer the same SSR for SNRs lower than 40 dB but in high 
SNRs, the RNN lA design provides a slight SSR improvement 
as compared with the NN lA design. The reason is that the 
later uses better surrogates of rank function than the former. In 
this system, all lA designs still outperform the conventional lA 
design, and the proposed lA schemes perform better than the 
WSLM lA algorithm. However, they perform worse than the 
ZFWS lA algorithms. The reason is probably that the perfect 
lA is not guaranteed although the considered system is proper. 



Fig. 2. The average SSR versus SNR for (18 X 12, 9, 3)® systems. 



Fig. 3. The average SSR versus SNR for (15 X 15, 9, 3)® systems. 

V. Conclusion 

In this work, we presented the lA schemes for a secure 
multiuser MIMO system in presence of an eavesdropper. 


We reformulate the design problem to minimize the rank of 
interference subspace and wiretapped signal matrices subject 
to the full rank of the desired signal matrices. To tackle with 
nonconvexity of rank functions, we introduce two surrogate 
functions, namely NN and RNN. Then, we developed iterative 
algorithms based on the coordinate descent approach to obtain 
suboptimal solutions of the precoding matrices and receive 
subspace matrices. Numerical results show that the proposed 
lA designs perform the same or better lA approaches based 
on minimizing power of interference and wiretapped signals in 
the system where the perfect lA can be achieved. In addition, 
both proposed lA algorithms outperform the conventional lA 
algorithm in terms of SSR. 

References 

[1] A. Mukherjee, S. A. Fakoorian, J. Huang, and A. L. Swindlehurst, 
“Principles of physical layer security in multiuser wireless networks: 
A survey,” IEEE Commun. Sun>eys & Tutorials, vol. 16, no. 3, pp. 
1550-1573, 2014. 

[2] M. Hanif, L.-N. Tran, M. Juntti, and S. Glisic, “On linear precoding 
strategies for secrecy rate maximization in multiuser multiantenna 
wireless networks,” IEEE Trans. Signal Process., vol. 62, no. 14, pp. 
3536-3551, Jul. 2014. 

[3] N. Yang, G. Geraci, J. Yuan, and R. Malaney, “Confidential broadcasting 
via linear precoding in non-homogeneous MIMO multiuser networks,” 
IEEE Trans. Commun,, vol. 62, no. 7, pp. 2515-2530, Jul. 2014. 

[4] G. Geraci, M. Egan, J. Yuan, A. Razi, and I. Codings, “Secrecy sum- 
rates for multi-user MIMO regularized channel inversion precoding,” 
IEEE Trans. Commun., vol. 60, no. 11, pp. 3472-3482, Nov. 2012. 

[5] H.-M. Wang, F. Liu, and X.-G. Xia, “Joint source-relay precoding 
and power allocation for secure amplify-and-forward MIMO relay 
networks,” IEEE Trans. Inf. Forens. Security, vol. 9, no. 8, pp. 1240- 
1250, Aug. 2014. 

[6] H. Jeon, S. McLaughlin, I.-M. Kim, and J. Ha, “Secure communications 
with untrusted secondary nodes in cognitive radio networks,” IEEE 
Trans. Wireless Commun., vol. 13, no. 4, pp. 1790-1805, Apr. 2014. 

[7] Y. Pei, Y.-C. Liang, L. Zhang, K. Teh, and K. H. Li, “Secure commu¬ 
nication over MISO cognitive radio channels,” IEEE Trans. Wireless 
Commun., vol. 9, no. 4, pp. 1494-1502, Apr. 2010. 

[8] O. Koyluoglu, H. El Gamal, L. Lai, and H. Poor, “Interference align¬ 
ment for secrecy,” IEEE Trans. Inform. Theory., vol. 57, no. 6, pp. 
3323-3332, Jun. 2011. 

[9] S. Sasaki, T. Shimizu, H. Iwai, and H. Sasaoka, “Secure communi¬ 
cations using interference alignment in MIMO interference channels,” 
in Proc. IEEE Int. Symp. Antenna. Propagat. (ISAP), Oct. 2012, pp. 
762-765. 

[10] T. T. Vu, H. H. Kha, T. Q. Duong, and N.-S. Vo, “On the interference 
alignment designs for secure multiuser MIMO systems,” submitted to 
lET Communications, 2015. 

[11] D. Papailiopoulos and A. Dimakis, “Interference alignment as a rank 
constrained rank minimization,” IEEE Trans. Signal Process., vol. 60, 
no. 8, pp. 4278^288, Aug. 2012. 

[12] H. Du, T. Ratnarajah, M. Sellathurai, and C. Papadias, “Reweighted 
nuclear norm approach for interference alignment,” IEEE Trans. Com¬ 
mun., vol. 61, no. 9, pp. 3754—3765, Sep. 2013. 

[13] S. W. Peters and R. W. Heath, “Interference alignment via alternating 
minimization,” in Proc. IEEE Int. Conf. Acoust. Speech. Signal Process. 
(ICASSP), Taipeo, Taiwan, Apr. 2009, pp. 2445-2448. 

[14] M. Fazel, H. Hindi, and S. Boyd, “Rank minimization and applications 
in system theory,” in Proc. IEEE Conf. Amer. Control (ACC), vol. 4, 
Boston, USA, Jul. 2004, pp. 3273-3278. 

[15] M. Grant and S. Boyd, “CVX: Matlab software for disciplined convex 
programming, version 2.1,” http://cvxr.com/cvx. Mar. 2014. 

[16] S. Bazzi, G. Died, and W. Utschick, “Interference alignment via min¬ 
imization projector distances of interfering subspaces,” in Proc. IEEE 
Int. Conf. Signal Process. Advances in Wireless Commun. (SPAWC), 
Cesme, Turkey, Jun. 2012, pp. 274—287. 



















